Do Rotations Beyond the Cosmological Horizon 
Affect the Local Inertial Frame? 



m 
o 
o 



> 
(N 

o 
m 
O 

cr 



X 



Jifi Bicak, 1,2,4 Donald Lynden-Bell, 2,1,4 and Joseph Katz 3,1,4 
1 Institute of Theoretical Physics, Charles University, V Holesovickdch 2, 180 00 Prague 8, Czech Republic 
^Institute of Astronomy, Madingley Road, Cambridge CBS OHA, United Kingdom 
3 The Racah Institute of Physics, Givat Ram, 91904 Jerusalem, Israel 
A Max-Planck Institute for Gravitational Physics (Albert Einstein Institute) 144-76 Golm, Germany. 

If perturbations beyond the horizon have the velocities prescribed everywhere then the dragging 
of inertial frames near the origin is suppressed by an exponential factor. However if perturbations 
are prescribed in terms of their angular momenta there is no such suppression. 

We resolve this paradox and in doing so give new explicit results on the dragging of inertial frames 
in closed, flat and open universe with and without a cosmological constant. 
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I. INTRODUCTION 

In a clearly written paper Schmid Q] describes how 
rotational perturbations of a spatially flat universe influ- 
ence the inertial frames. He finds that such influences are 
attenuated by an exponential Yukawa factor whenever 
the perturbation lies beyond a 'horizon'. He expressed 
his results in terms of a quantity that he calls the energy 
current J e . What corresponds to Schmid's energy current 
J e in our calculation is a quantity J s with the dimensions 
of angular momentum. However his result appears to dis- 
agree with our earlier result |2f (hereafter LKB) that in 
a flat universe the rotation of the inertial frame, u>, due 
to any system of spheres with small rotations about a 
center is given by 



w(r,t) = -j- 



^J(t,<r) 



1 dJ 
— o —r— dr 

r' 3 dr' 



(1.1) 



where J(i, < r) is the total angular momentum within 
the sphere of proper radius 3] r — a{t)r. This expres- 
sion demonstrates how angular momenta at all distances 
contribute and shows no exponential cut-off and no in- 
fluence of any horizon. 

Both results agree, however, that inertial 'influences' 
may be expressed instantaneously i.e. with no light 
travel-time delay. This is because they follow from the 
constraint equations of General Relativity with an appro- 
priate mapping onto an unperturbed universe to provide 
a suitable gauge. 

The results are in apparent contradiction. However, 
when the details of both calculations are examined it is 
evident that the contradiction lies in the attribution of 
different causes for the effect. Schmid's 'energy current' 
is considered by him as the source of the rotational drag- 
ging of inertial frames lu. Schmid's quantity does not 
obey a conservation law but for a rotating sphere it can 
be directly expressed in terms of its angular velocity f2, 



density p, pressure p and proper radius r as 

J t = 2ir (p + p)nsm 3 e r A d9dr, (1.2) 



fl is the 'coordinate' angular velocity, not that measured 
relative to the local inertial frame which, as we show be- 
low, is the quantity needed in the angular momentum 
conservation law. The contribution to the conserved an- 
gular momentum is 



2ir(p +p)(Cl - lu) sin 3 Or 4 d9d 



(1.3) 



The difference comes about mathematically because the 
perturbed metric is not diagonal but in spherical polar 
coordinates is 



ds 2 = dt 2 



-a 2 (t){ 



dr 2 



dO 2 



-sin 



'0{dLp 2 -2LjdLpdt) j, 

(1.4) 

where Lu(t, r) gives the small rotation of the inertial frame 
due to the rotational perturbations everywhere. 

For a perturbed 3-flat universe, angular momentum 
conservation is given by the equation 



d 



(1.5) 



where rf is the angular Killing vector of the background 
(flat) space corresponding to the particular component 
of angular momentum considered. 

Thus the conserved quantity is (the minus sign comes 
from the signature we use in the metric) 



J = 



_rpO 



-gd 3 x. 



(1.6) 



Since the ip component of T^ v is brought down in this 
expression, it is not merely the motion of the fluid that 



is involved in T°„ via its contribution 



but also the 



off-diagonal metric component go v which depends on u) 
at the position of the source [see our metric l|1.4Jl ] ■ As we 
see from equation (|1.1|) , we regard the conserved angular 
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momentum J as the source of the dragging of inertial 
frames, and this was the quantity we used in LKB. 

Schmid's work for a spatially flat universe is more 
general than our work published so far, since he con- 
siders all vector perturbations, nevertheless we treated 
closed and open universes as well as flat ones and in- 
deed from a Machian viewpoint it is the closed universes 
that are more interesting by far. We also considered 
all spherical but inhomogcncous Lemaitre-Tolman-Bondi 
universes with rotational perturbations that were con- 
stant on spheres. Finally we looked into the problem of 
the rotation of inertial frames induced by spheres of given 
angular velocities, rather than given angular momentum. 
This is a special case of Schmid's problem but general- 
ized to closed and open universes. In our discussion we 
wrote down the equations governing u>(t, r) when Cl(t, r) 
was given and showed how they could be solved. We car- 
ried out the detailed solution only for the static closed 
Einstein universe (LKB Appendix A). 

Schmid's beautiful result that the dragging is exponen- 
tially suppressed when a sphere of given angular velocity 
is outside his horizon, has stimulated us to work out all 
our solutions in detail for all FRW universes. Barring fac- 
tors of a(t) that Schmid seems to have omitted in error, 
we fully confirm his result for a flat universe. Thus we 
have the fascinating paradox that while spheres of given- 
angular velocity have their dragging exponentially sup- 
pressed if they are outside the horizon, nevertheless the 
dragging of spheres of given angular momentum suffers 
no such suppression! How can this be! 

In the prescribed angular momentum problem one may 
consider (for an open or 3-flat universe) having only one 
spherical shell of finite thickness having angular momen- 
tum. The gravity of this source will induce a rotation of 
inertial frames everywhere. The fluid at all other places 
will respond inertially and start to rotate so that Q = lu 
everywhere except on the original shell. Thus in the pre- 
scribed angular momentum problem we give one shell 
angular momentum, sit back and watch. We see how 
the inertial frames are affected everywhere else merely 
by watching the rotations of all other spheres. 

The prescribed angular velocity problem needs more 
organization in the creation of the initial state. If we 
start one thick shell rotating at the prescribed rate then 
all the other will start moving so as to keep up with the 
induced rotation of the inertial frames. If the prescription 
is to have just the one thick shell rotating and none of the 
others we shall have to stop them. In doing so we have 
to give them negative angular momentum to keep f2 zero 
even though the inertial frame is rotating at to. When 
in Schmid's problem the perturbation in angular veloc- 
ity is considered as confined beyond his horizon he shows 
that the rotation of inertial frames is exponentially sup- 
pressed near the origin. The prime reason is that in order 
to keep the motion confined, the intervening spheres have 
to be given backward angular momentum to stop them 
from following their inertial tendency of rotating at the 
inertial rate lu. The influence of all the backward angular 



momentum (of non-rotating spheres!) rather effectively 
cancels most of the rotation of the inertial frames induced 
by the original shell. Hence the suppression of the effect 
is due to all the negative angular momentum that was 
supplied to keep the other spheres from rotating! The re- 
maining suppression is due to the rotation of the inertial 
frame at the original shell itself: lu there is a fraction of 
£1 so that — lu, on which the source depends, is less. 

There is a long history of treating dragging effects 
within spheres starting with Einstein's treatment using 
an early version of his gravitational theories. Within 
General Relativity the early works of Thirring £| and 
Lense and Thirring 0] were later generalized to deep po- 
tential wells by Brill and Cohen This raised ques- 
tions as to whether the dragging would be perfect within 
a black-hole's horizon. We believe that the first paper to 
remark on the apparent instantaneity if inertial frames 
is the pioneering paper of Lindblom and Brill Q on in- 
ertia in a sphere that falls through its horizon. More re- 
cently we explored observational effects seen within such 
a sphere Q and gave an example of strong linear drag- 
ging in a rapidly accelerated charged sphere 0. Strong 
cosmological perturbations in a weakly rotating sphere 
surrounding a void were treated by Klein [Io| and in 
greater detail by Dolezel, Bicak, and Deruelle [llj who 
also discussed how an observer within such a cosmologi- 
cal shell views the world outside. 

We owe a debt to Schmid as his work stimulated us to 
work out the consequences of our solutions Q in much 
greater detail and, without that, we would never have 
raised, let alone understood, the delightful paradox em- 
phasized above. In particular we have now investigated 
thoroughly the problem when lu is to be solved for with 
the angular velocities given everywhere at one cosmic 
time. Previously we had concentrated on the problem 
with the angular momenta given. While both are impor- 
tant problems we strongly believe that it is the latter that 
is of dynamical importance in formulating Mach's Prin- 
ciple. It can nevertheless be argued that the apparent 
agreement between the angular positions of quasars at 
different epochs and the inertial frame defined by using 
the solar system as a giant gyroscope stimulates Machian 
ideas. While it is the angular momentum that is impor- 
tant for the physics it is the apparent kinematical agree- 
ment between the angular velocity of the sphere and the 
inertial frame that is observed. In this sense the prob- 
lem with given angular velocities may be closer to Mach's 
original and it is unclear how distant observations could 
measure the true angular momentum of a sphere includ- 
ing its dragging term, while its angular velocity is more 
directly observable. However, see [ill f° r the complica- 
tions of light bending. Unfortunately the problem of the 
observed agreement of frames is not that either Schmid 
or we have addressed since both our treatments relate in- 
stantaneous quantities at the same cosmic time whereas 
observers use no such world map (except in the solar 
system) but a world picture in which distant objects are 
seen as they were long ago on the backward light cone. 
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It seems unlikely that an exact causal relationship exists 
between proper motions of masses on our past light cone 
and our local inertial frame, since, at any cosmic time 
the inertial frame 's rotation has contributions from ob- 
jects that were never in our past light cone. Of course 
such objects will no doubt have been seen by some alien 
and the Copernican principle would suggest that the ap- 
parent agreement of the kinematic and inertial frames 
here will be repeated there. What is under discussion 
above is the influence of distant bodies on the local in- 
ertial frame. This is quite distinct from a comparison 
of the dynamics of the solar system with its kinematics 
relative to distant quasars (as seen on a hundred years of 
past light cones), from which the rotation of the inertial 
frame is computed. 

Beside the resolution of the apparent contradiction 
with Schmid the main contributions of this paper are 
the following. 

Section 2. The derivation of the equations govern- 
ing general perturbations and a brief introduction to 
Machian gauge conditions which allow the separation of 
the (hok) vector perturbations equations from the others. 
The discussion of the equations of motion that must be 
obeyed if the contracted Bianchi identities are to be sat- 
isfied. As a consequence when axial symmetry is imposed 
each ring of fluid preserves its angular momentum. This 
section concludes with basic equations for odd parity axi- 
ally symmetrical perturbations from which the remainder 
of the paper is derived. 

Section 3 derives the explicit expressions for rotation 
of inertial frames in terms of the angular momentum dis- 
tribution at any one time. This is done for all FRW uni- 
verses with k = ±1 or but the simplest case is solved 
in this section with f2 constant on spheres at the time 
considered. This corresponds to odd-parity vector I = 1 
perturbations with Q independent of 9. In the following 
paper 0] (Paper II) we allow for general dependence. 
With the integrals evaluated at fixed cosmic time and 
with the constants c and G restored we have the follow- 
ing results for Q(r) at fixed time (for the derivation of 
the vector forms below see 0): 

For k = 0, r = a(t)r, 



For k 



2G 



J(< r)i 



°° dJ -3 j 
— r 3 dr 

dr 



(1.7) 



Notice that lu cx 
For k = 1, 



[a(t)] 3 since J is conserved. 



2G 



J(< x)w( x ) 



dJ 
dx 



■u {t), 



(1.8) 

and W( X ) = 



here LOa(t) is undetermined, r 

cot 3 x + 3cotx- The arbitrariness of C3o(t) is intimately 
connected with Mach's principle. The physical J involves 
(51 — uj) and does not change for rotating axes as it in- 
volves a difference, see Q and below. 



2G 



J(< x)w(x) 



f°° r] J 



(1.9) 



where W(x) = coth 3 x~ 3cothx + 2, and W has an extra 
2 so it tends to zero at x ~~ * 00 thus ensuring that the 
boundary condition uj — > is obeyed. When contribu- 
tions from a 6 dependence of f2 are included these results 
are supplemented by 9 dependent terms that average to 
zero on spheres. More general results are given in the 
accompanying Paper II. 

Section 4 gives explicit solutions for the rotations of in- 
ertial frames for the same special forms of perturbations 
as in Section 3 but now it is the angular velocities of the 
different spheres that are given rather than their angular 
momenta (this is closer to what might be observed but 
cf. earlier discussion). We define A by 



A 



2 K a 2 (p + p) =A(k-a 2 H), 



(1.10) 



8ttG/c 4 



Sit in geometrical units used in the 



following, the dot denotes d/dt, H = a/ a is the Hubble 
constant. The second relation in l|1.10(l follows from the 
combination of the background Einstein's equations for 
any p, p, k and also for any value of the cosmological con- 
stant A. The rotation of inertial frames near the origin 
due to an fl distribution at large z' = Xr is for k = 



ui(r) 



z' 2 e~ z n(z')dz', (1.11) 



which shows Schmid's exponential attenuation e z . At 
the perturbation itself, close to zq, we find for z' large: 

For k = 1 we give the results near the origin and at 
the perturbation when A 2 > 4. When A 2 < 4, which can 
occur when a A-term is present, there is no exponential 
in the expression. It is assumed that exp(vA 2 — 4%) is 
large at the source. With r = sinx 



1 ( A 2 x 2s 

^ =3 1 + -nr 



(1.13) 



AVA 2 -4 er^^x' sin 2 ( X 'Mx') dx'- 



We have assumed exp(vA 2 — 4 x) 3> 1 for x — tt, x' , and 
7r — x' ■ At the 'source' 



^(Xo) 



I A 2 -3 / sm\- 



-v / A^4| x '- X0 | 



2 Jo \/A 2 - 4 V sin Xo, 
Similarly for k = — 1: r = sinhx, 

-(x)=K i+ ¥' x 



(1.14) 
"(xW- 



(1.15) 



AVA 2 + 4 e-^^x' sinh 2 (x'Mx')dx', 
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and at the source 

w(xo) 
1 



/ sinh x' 
2 Jo VA 2 +4 V sinh Xo 



,-VX r +S\x'-xo\ 



(1.16) 



We emphasize that all of the above relationships are 
true at any given instant, but that both the angular mo- 
mentum distribution and the angular velocity distribu- 
tion at later instants are related to those at earlier times, 
so can not be given independently of those given at an 
earlier epoch. In axial symmetry the angular momentum 
distribution follows the motion of the perfect fluid but, as 
the angular momentum is first order and the movement 
across the background is of first order, the product can 
be neglected. Thus to first order the angular momentum 
density can be considered as painted on the background. 
This is not true of J s which is not conserved and nor is 
it true of the angular velocity O. In both cases to find 
the time evolution one must appeal to the equations of 
motion which, in axial symmetry, leads back to local con- 
servation of angular momentum density. Only by use of 
its conservation can one find how Q and J s can evolve 
consistently with Einstein equations (i.e. with the con- 
tracted Bianchi identities). In this sense the given angu- 
lar momentum problem is far more physical than either 
Schmid's problem or the given f2 problem to which it is 
equivalent. The time evolution of u> and Q are derived 
and discussed in Section 5. 

In a paper that has long been in gestation we give 
a discussion of those gauges in which the Machian rela- 
tions of the local inertial frames to the motions of distant 
masses can be expressed instantaneously at constant cos- 
mic time. In that paper we derive all equations that gov- 
ern all perturbations. All can be solved using harmonics 
in the 3-space of constant time. However harmonics are 
not as informative as Green's functions so in the follow- 
ing paper we integrate the relationships between the 
rotations of the inertial frames and the angular momen- 
tum density for all axially symmetrical odd-parity vector 
perturbations, called usually "toroidal" perturbations in 
astrophysical and geophysical literature. These results 
allow — u>, which enter the angular momentum den- 
sity, to be any function of (r, 6) but independent of (p. 
However, since the background is spherically symmet- 
ric, non-axisymmetric perturbations can be generated by 
re-expanding axisymmetric perturbations around a new 
axis, and taking the component with the new e lTav as the 
component with general m. 



where the background metric is used to move indices 
and the time-independent part of the spatial background 
metric [i, j,k = 1,2, 3] is used to define the 3-covariant 
derivative V k and V fe = f kl Vi- 

In one of the standard coordinate systems the back- 
ground FRW metric reads 



ds 2 = dt 2 



dr 2 



1 — kr 2 



+ r 2 (d6> 2 +sin 2 9dip 2 ) 



(2.2) 



where in positive curvature (closed) universe (k = +1) 
r £ (0, 1), in flat (k = 0) and negative curvature [k = — 1) 
open universes r G (0,oo), and 9 e (0, if), <p S (0,2tt). 
We shall also employ hyperspherical coordinates 

ds 2 = dt 2 - a 2 [d X 2 + r 2 {d6 2 + sin 2 6>d<^ 2 )] , (2.3) 

with r = sinx, X: sinhx for k = 1, 0, — 1. 

In a completely general gauge for general perturbations 
h^v, the (momentum) constraint equation, SGl = kST®, 
turns out to be 



V 2 h 



Ok 



kh ok 

z o 

= a 2 nST°, 

where the dot denotes d/dt 



o 2 V fc Vj/ij 



-aVt/C 



a T fc 



K 



■^Hh 00 



is the perturbed mean external curvature of t 
slices, H — a/ a is the Hubble constant, 



% = -V, 



(2.4) 



(2.5) 



constant 



(2.6) 



Notice that equation (|2.4() is independent of the choice 
of the cosmological constant A because we perturbed 
"mixed" components of G°. Other perturbed Einstein's 
equations will not be needed here. Since, however, we 
are interested primarily in perfect fluid perturbations we 
shall also consider the perturbed fluid equations of mo- 
tion, i.e. the perturbed Bianchi identities 



-Hhoo 



(8p)' + 3H(5 P + 5 P ) + 
( P + p)V k (h k Q +V k ) + (p + p) 

and 

4 [a 3 ( P + p)(a 2 f km V m -ho k )] + 



a 



(2.7) 

= 0, 

(2.8) 



II. THE EQUATIONS TO BE SOLVED 

We write the perturbed FRW metric in the form 

d« 2 = {g^ + h^&xW (2.1) 
= dt 2 - ^(tjfcjdx'dx 1 + h^dx^<dx v , 



Vfc<5p + -(p + p)V fc /i o 



0, 

where V k = ^ ~ SU k and V k = -a 2 f k jV j is 
the fluid (small) velocity. The perturbed fluid energy- 
momentum tensor components entering the constraint 
equations (|2.4|l read 



ST° = (p+p)(h ok +V k ) 



(p+p)(h 0k -a 2 f km V m ). (2.9) 
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There have been various choices of gauges used in the 
literature, in particular the synchronous gauge (/loo = 
hok = 0). In order to understand the effect of drag- 
ging of inertial frames, in particular its 'instantaneous' 
character, it is convenient to use gauges — we call them 
'Machian' — in which the constraint equations, and still 
another (combination of) the perturbed field equations 
are explicitly the elliptic equations. In order to achieve 
this it is first useful to choose coordinates on t — constant 
slices such that the spatial harmonic gauge conditions are 
satisfied, i.e. T k = 0, where T k is given in (|2.6[l (in nu- 
merical relativity T k = is frequently called the 'minimal 
distortion' shift vector gauge condition). Next, it is con- 
venient to choose the time slices so that, for example, the 
perturbation of their external curvature vanishes: K, = 0, 
K, given by (|2.5[) (so called 'uniform Hubble expansion 
gauge'). Under these gauge conditions (which determine 
the coordinates in a substantially more restrictive way 
than e.g. the synchronous gauge) the constraint field 
equations l|2.4|l become the elliptic equations for just the 
components h ok , no other huv enter. 

Until now we considered general perturbations in the 
chosen gauge. Hereafter, we assume the vectors hot, V k 
to be transverse, 

V% fe - 0, V k V k = 0, (2.10) 

so that also V k 8T^ = 0. If 12.11)11 is not satisfied, we 
can apply V fc to equation (|2.4I) . find the elliptic equation 
for the scalar V k tiQ k , solve it and substitute back into 
(|2.4|) where the third term on the left hand side could be 
viewed as the source together with ST®. Since, however, 
the longitudinal parts do not contribute to the dragging 
of inertial frames, we assume equations l|2.10|l to be sat- 
isfied. 

The constraint field equations l|2.4|l with our choice of 
gauge JC = T k = [cf. equations l|2.5|l and l|2.6|l ] thus 
become 

V 2 h 0k +2kh 0k = 2a 2 K 8T® = 2a 2 K(p+p)(h 0k -a 2 f krn V m ), 

(2-11) 

where for the perfect fluid ST® is substituted from equa- 
tion (|2.9|) . This is our basic equation to be solved at a 
given time t = constant, with either ST® or V k given. 
The Bianchi identities (fluid equations of motion) de- 
termine the time evolution of perturbations, the scalar 
equation (|2.8|l for 5p, whereas the vector equation (|2.9() 
governs the evolution of the term 

Jk ee a 3 (p +p)(a 2 f km V m - h Qk ) = -a 3 ST®. (2.12) 

In the following we shall often express the background 
time dependent term a 2 (p + p) by using equation Hl.lOfl • 
Consider first the flat universe (k = 0). In Cartesian 
coordinates x k used by Schmid pj , the 3-metric f k i = 5 k i , 
and p.llfl becomes 

V 2 h 0k = 2a 2 k6T% = 2a 2 n(p + p)(h ok + V k ), (2.13) 

where V 2 is the flat-space Laplacian. Substituting from 
equation (|1.1U|) with k = in the first term in the r.h.s. 



of equation l|2.13|l . we get 

V 2 /i fc = -4a 2 Hh 0k - 2Ka i (p+p)V k . 



(2.14) 



Now comparing our general form of the perturbed 
FRW metric with the perturbed metric (5) in Schmid's 
work (and bewaring of the opposite signature), we see 
that ho k = — a/3fc(Schmid). Considering (p + p)V k (de- 
noted by J e in Schmid) as the source, the equation (|2.14|) 
written for Schmid's (3 k becomes 



-V 2 /3 fc - 4a 2 H/3 k = -2 K a 3 (p+p)V k . 



(2.15) 



This is Schmid's basic equation (14), up to the factors a 2 
and a 3 which in Schmid's equation (14) are missing but 
this does not change significantly Schmid's conclusions. 

When ST® is given, the solution of equation (|2.13|l is 
given as the Poisson integral over the source. If, how- 
ever, the matter current is given, equation l|2.14|l can be 
written as 

V 2 h 0k - X 2 (t)h Qk = -2na 4 {p + p)V k , (2.16) 
with (k = 0) 

A 2 = -4a 2 iJ. (2.17) 

Usually (e.g. in the standard Friedmann models) H < 0, 
so A is real. The three equations l|2.16|l are, as empha- 
sized by Schmid, of the Yukawa-type. The Green's func- 
tions are given by 



G(x,x') 



(2.18) 



4-7T \x — x'\ 

the well-behaved solution of equation (|2.18(l is thus 

hok = -^™ 4 (P + P)J W) 6 ,. da/. (2.19) 

Clearly if the perturbation V k (x') is located at \x — 
x'l^A" 1 = l/2a\/^H, i.e. beyond the L H radius' R k = 
2(— H)~2 in Schmid's terminology, the vector ho k which 
determines the dragging of inertial frames is exponen- 
tially suppressed around the origin. Although we thus 
verified the interesting conclusion of Schmid, we do not 
resonate with his view that " because of the exponential 
cut-off... there is no need to impose 'appropriate bound- 
ary conditions of some kind'....". The Green's function 
in (I2.18|l with the '+' sign in the exponential is also the 
solution of equation l|2.16|l with a (5-function source but 
one discards it by demanding 'reasonable' boundary con- 
ditions at infinity. 

From the Machian viewpoint the closed universes are 
of course preferable. There is, however, no vector Green's 
function available for equation (|2.11|l with either <5T" or 
V k considered as a source. In order to understand how 
Schmid's conclusions get modified in curved universes 
and to generalize our previous work [2j which analyzed 
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perturbations corresponding to rigid rotating spheres in 
the FRW universes, we shall study all axisymmetric, odd- 
parity Z-pole perturbations corresponding to differentially 
rotating 'spheres'. We now derive the basic equations for 
such "toroidal" perturbations. Their solutions, in par- 
ticular for I > 2 and closed universes, require special 
treatment. These solutions are analyzed in the following 
Paper II. 

In spherical coordinates [as in the FRW metrics l|2. 2J1 . 
(|2.2|) ]. the only non- vanishing vector components are 
hoip(t,r,9) and V v (t,r, 9). [For the general axisymmet- 
ric even-parity vector fields V v = 0, whereas V r (t, r, 9) 
and Vg(t, r, 9) are non- vanishing, the same being true for 
honhoe]- There is now just one non-trivial constraint 
equation in (|2.11() to be satisfied: 



Otp 



2kh 0v = 2a 2 n5T°, 



(2.20) 



in which V = / VfcVj, with / being the inverse to fki 
given by FRW metric l|2.2|l (recall - see (I2.2|l - that f k i is 
positive definite, without factor a 2 ). Calculating V 2 /io v 
explicitly, we find equation H2.20J1 to take the form 



(1 - kr 2 ) 
„ d 



1 

sin ( 



d 2 

dr 2 


kr ir 

or 


( 1 




Vsinf 


189 



h 0<p + (2.21) 
io<p + 4:kh 0ip = 2a 2 nST°. 



Before solving this constraint equation it is interest- 
ing to notice what the perturbed equations of motion 
(Bianchi identities) say for axisymmetric odd-parity per- 
turbations. Equation l|2.8Jl in our gauge choice (with 
K, = 0) and transverse character of hok,Vk is a simple 
evolution equation for Sp. The vector equation 1|2.9|) for 



indices 1, 2 (a; 1 



r, x 



9) turns into the well known rel- 



ativistic equilibrium conditions for perfect fluids, Vfc<5p = 
— (p+p)Vk(hh-oo) (see e.g. ^3|). In the following the cru- 



cial role is played by equation l|2.9|l for index 3 
Since in axisymmetric case V v Sp = 0, V v /ioo = 
comes 



[a 3 (p + p) (a 2 r 2 sin 2 0V V - h 0v )] 



[a 3 ST^]' = 0. 



o. 



x = ip). 
0, it be- 



(2.22) 



(2.23) 



This is the conservation of angular momentum of each 
element of each axially symmetrical ring of fluid. The 
total angular momentum in a spherical layer (xi,X2) is 
given by 

J(X1,X2)=- [ X d X 1*69 / d<p^STy, (2.24) 
Jxi Jo Jo 

where r/^ = (0,0,0,1) is the rotational Killing vec- 
tor, the background metric determinant g = g^ = 
—a 6 r 4 sin 2 9, r = sinx,Xi sinhx for respectively k 



■1.0. 1 as in equation 12.3fl . Integrating over ip we 
have 

J(Xi,X2) = -2tt d X I d9a 3 r 2 sm96T° 
J X i Jo 

= 2tt / X2 d x f d6j(6, X ,t), (2.25) 
J X i Jo 

where j(6, x,t) is the (coordinate) angular momentum 
density. Hence, the Bianchi identity (|2.23|) can be written 

as 



(2.26) 



This is important for studying the time evolution of the 
hok and Vk perturbations. 

Defining the fluid angular velocity 



dt 



(2.27) 



we get 



a 2 f vv V ip = ~a 2 r 2 sin 2 9 fi(i, r, 9). (2.28) 



Writing similarly 

h 0v = a 2 r 2 sin 2 9 u(t,r, 9), (2.29) 

the only non-vanishing component of ST® becomes 

<5T° = {p + p)a 2 r 2 sin 2 9(lu - fi). (2.30) 

The angular momentum density conservation law l|2.23|l , 
resp. <|2.26[) . turns then into the simple evolution equa- 
tion 



[a 5 (p + p){cj -fi)]" = 0. 



(2.31) 



Let us now return back to the constraint equation 
(|2.22() . The second term on its left hand side suggests 
the decomposition into the vector spherical harmonics. 
It should be emphasized that, in contrast to standard 
practice in the cosmological perturbation theory where 
perturbations are decomposed into harmonics in all three 
spatial dimensions (see e.g. [13), we decompose in the 
usual coordinates 9,p on spheres only, and assume axial 
symmetry (spherical functions Yj m having m = 0). Thus, 
we write (Yi ,e = dgY m ) 



r 2 ^^(t,r)sin^ , e , (2.32) 
i=i 

oo 

V v = -aV^Oi^Osin^, (2.33) 



i=i 



and 



5T° = a 2 (p +p)r 2 J2(ui - Q t ) sine Y l0>e 
i=i 

oo 

= XW^W^o.e- (2-34) 
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Substituting these expansions into equation (|2.22|l and 
using the orthogonality of functions sin# Y;o,0 for differ- 
ent Z's, we obtain the 'radial' equation for each I: 



d 2 



d 



(r 2 LUi) - 1(1 + 1)loi + 4kr 2 uji 



dr 2 dr 

2a 2 n(p + p)r 2 (uj l - Hi) = \ 2 r 2 {iJi - Qi), (2.35) 



where we used equation l|1.10|l . It is easy to convert the 
last equation into the form 



d 



r z or 



VI - kr 2 — (r 2 w ; ) 
or 



l l±H LJl -Aku Jl = X 2 (n l -Lu l ). (2.36) 

For I = 1 (and the background pressure p = 0) this 
equation coincides exactly with equation (4.32) in LKB. 
In the language of the present paper, in LKB we ana- 
lyzed dipole (1 = 1) axisymmetric odd-parity perturba- 
tions. With I = 1, Y W g = — A/3/47T sin (9, so that putting 
lj = — y3 /47r Q = —y/3/4:TT Hi— i, we recover 

h 0v = a 2 r 2 sin 2 9 u(t,r), V v = -a 2 r 2 sin 2 9 Q(t, r), 

(2.37) 



5T° = a 2 (p + p)r 2 sin 2 9 (w~tt), 



which corresponds to the rigidly rotating spheres in the 
FRW universes considered in Section 4.4 in LKB, and, 
for n(t,r) given, analyzed in detail in Section 4 in the 
following. 

Consider first the case k = 0. Equation H2.36|l can be 
written with the angular momentum density (5T®)i as a 
source, 

1 dftduA 1(1 + 1) - 2 2 2/ W<h 

(2.38) 

If the fluid angular velocity is taken as a source, the equa- 
tion reads 



r 4 dr 



dr J 



-A O;, 



(2.39) 

where A 2 = —Aa 2 H = 2na 2 (p + p) by using equation 
(QUI) with k = 0. 

In the case of spatially curved (k ^ 0) backgrounds it is 
advantageous to write r 2 = k(l — p 2 ), i.e. p = y/1 — kr 2 
to obtain 



1 



d 



[k(l- /i 2 )] 3 / 2 <fy 



{[*( 



l-/i 2 



is/agj^il ?(/ + !) -2 
J dp J A:(l-M 2 ) 



2/v 



A;(l-/i 2 ) 



(5T°)z. (2.40) 



The substitution 



w, = [k(l - p 2 )}-^ 4 



(2.41) 



turns equation (|2.40|) into the Legendre equation for u>i 
with (5T®)i as the source: 



d_ 

dp 



2\ au} l 



¥2 +1 
2k 



(l + \) 2 



to i 



k(l ~ p 2 ) 

w^w (KV (2 - 42) 



Finally, considering the fluid angular velocity as the 
source, we can write the last equation again as the Leg- 
endre equation with a more complicated degree: 



d_ 

dp 



2 , duji 



kv(v + 1) 



(I 



1)2 
2> 



LJl 



k(l-p 2 )_ 

= -K l = -\ 2 n l [k(l-p 2 )f'\ (2.43) 



where 



1 



v + - = 4-2kKa 2 (p+p) = 4-kX 2 = 4ka 2 H. (2.44) 



The degree v of the Legendre equation does not depend 
on I. For I = L equation (|2.43l) goes over into equation 
4.35 in LKB 



III. SOLUTIONS FOR w WITH GIVEN 
ANGULAR MOMENTUM DISTRIBUTION 

We shall start by making more explicit the solutions 
obtained in LKB which are the I = 1 odd-parity vector 
solutions of the general problem. In such modes each 
sphere rotates with no shear but it expands (or contracts) 
with the background and as it does so its angular velocity 
changes (see Section 5). 

k = 

The equation to be solved is (|2.38|l with 1 = 1, this is 
4.33 LKB 



1 d 



i duj s 



multiplying up by r 4 this takes the form 







i duj 



6 dJ(< r) 



so 



dr V dr 

duj 

dr 



dr 



6 J 



(3.1) 



(3.2) 



(3.3) 



the constant of integration is zero since J(< r) is zero at 
r = where du/dr must vanish. Integrating again and 
insisting that lu — > at oo we find 



— A Ar=2a 



J(<r) 



°°dJ ,_ 3 , 
-— r dr 
dr 



2-rrr' 2 sm9(-ST°„)d9dr' 



'0 jo 



J 2Trr'~ 1 sm9(-5T°)d9dr', 



(3.4) 
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where we have used (|2.25|) to define J(< r) in terms of 

k = i 



The equation to be solved is (|2.4UI) with I = 1 which is 
4.34 LKB [3 



so that 



, duo ' 



6 dJ 



i duo 



GJ 



dfi a 3 (l — /j, 2 ) 2 



(3.5) 



As before there is no integration constant for the same 
reason. We now write [i = cosx, then \ is the normal 
cosmic radial angle and 



6 J 



da; 



9% a 3 sin 4 x 



(3.6) 



Now 



dx 
sin 4 x 



(cot 3 x + 3cot x ) =--W{x). (3.7) 



Hence 



uj = 2a 



WJi< K )+ j w^d x ' 



wo, 



(3.8) 



where 



,/ = 



27ra 3 [r(x')] 2 sin 2 8(-6T°) d9d X '. (3.9) 



o Jo 



Just as in the last case W diverges at x = like x~ 3 , 
however, the angular momentum of spheres near the ori- 
gin is sufficiently small to make the WJ tend to a con- 
stant as x tends to zero. It is shown in LKB that the 
condition of convergence of the second integral at x — ^ 
is that the total angular momentum of the universe is 
zero. If that condition is fulfilled and Q — lu is regu- 
lar near x = n then the integral converges. If the total 
angular momentum is not zero then the integral for lo 
diverges at x = t- Thus for lo to be finite at x = ^ 
the total angular momentum must be zero in the closed 
universe. There is no way of fixing u)q because there is 
no standard of zero rotation, as there is for the infinite 
universes. Indeed, according to Mach a description of 
the world in rotating axes is just as good in principle as 
a description in non-rotating ones. Note that the source 
f2 — lo does not change when the axes are rotating since 

and lo acquire the same constant lOq. An absolute ro- 
tation can arise only from spatial boundary conditions 
which do not occur for closed universes. 



The equation to be solved is l|2.40jl with 1 = 1. Multi- 
plying through by (fi 2 — l) 3 / 2 we obtain 



d_ 



1)5/2^1 
' dfl! 



6 dJ 
a 3 d/x ' 
(3.10) 



so on integration and division 

Writing /i = coshx to introduce the natural radial vari- 
able of hyperbolic space, this becomes 



8lo 
dx 



6 J 
a 3 sinh 4 x 



(3.12) 



Integrating again and insisting that lo — > at infinity we 
use the integral 



v 7 = --(coth 3 x-3coth X + 2)^--iy(x), 



sinh 4 x! ^ 
and on integrating by parts we obtain 



lo = 2a 



ir./K-\)+/ w(x')^d x ' 



(3.13) 



(3.14) 



where J is the same as in (|3.9|l with r = sinhx- We 
have chosen the above definition of W so that W — > 
at infinity; so no constant of integration is needed to 
incorporate the boundary condition that lo — > 0. 



IV. SOLUTIONS FOR lo WITH GIVEN Q 

The method of solution was outlined in LKB but here 
we work through all the details starting with the simplest 
case. 

k = 



The relevant equation to be solved is H2.39JI with 1 = 1, 
equation 4.33 in LKB, rewritten as 



_L__5 

r 4 dr 



{ du> 
Or 



X 2 lo = -\ 2 n. 



(4.1) 



Here A 2 = 2a 2 K,(p+p) > 0. A _1 a has the units of a length 
and we shall call it, following Schmid 0, the distance 
to the horizon. In dimcnsionlcss comoving coordinates 
this corresponds to r = A -1 . We write z = Xr and 
doojdz = lo' . Then equation (|4.1() reduces to 



— LO = — Q. 



(4.2) 



The corresponding homogeneous equation is Bessel's 
equation for z -3 / 2 J 3 / 2 (iz), which has real solutions lo = 
X and lo = K, where X = z~ 3 / 2 I 3 / 2 (z) and K. = 
z- 3 / 2 K 3/2 (z). For small z, X -► + z 2 /10); 

JC — * \/tt/2 z~ 3 . For large z, X 




We use the method of variation of parameters to solve 
the inhomogeneous equation with boundary conditions 
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that u> tends to zero at infinity and to a constant at the 
origin. We thus obtain 



To keep P 3 ^ 2 (cosx) and Qt /z (cosx) rea h we use (1 — 
/i 2 ) _ 2 = (sinx) 1 in place of (/i 2 — 1)" in the recurrence 
relation 8.5.1 of Abramowitz and Stegun 01 (this merely 
multiplies the results by —i). 



3/2, 



K(z) / (z') 4 l(z')n(z')dz' 



For the solutions near the origin with sources that are 
not so close, we may neglect the first term and then for 
small z, 



{z'fK,{z')n{z')dz'. (4.4) 



When the source f2 is beyond the horizon z = 1, i.e. 
z' > 1, 

= ~ (l + Q J™{z'fe- Z 'n{z') dz'; (4.5) 

so for a source localized in ro(l ± A) with A <C 1/A, 

= J (l + ^) (Ar ) 3 e- Ar ^2A, (4.6) 

which clearly shows the exponential decline of influence 
remarked on by Schmid 1] . When Q is concentrated near 
Zq, in zq ± AA, then with zq 3> 1 and f2 = f2 we get 



w(zo) 



2 Jo 



oo / / \ 2 



zq 



-\ z '- Zc >\VL{z')dz' ~ AAfl (4.7) 



Thus a< the source the inertial frame rotates at AAf2 and 

U-u = (i- aa)H. 

We now turn to the solutions for a closed universe. 
fc = 1 

The relevant equation is Legendre's equation for ZJ = 
(1 — /i 2 ) 3 / 4 w with an inhomogeneous term written below. 
This is LBK equation 4.35 and the same as the equation 
of Section 2 of this paper specialized for I = 1: 

o r,„ o,^u>^ ( , ^ (|) 2 i_ 



(4.8) 



(l-/r 

# = A 2 r>(i- M 2 ) 3 / 4 , 



where (i/ + i) 2 = 4 - A 2 as in (|2.44() . Since k = +1 the 
space is hyperspherical and the convention is to write 
fi = cos\ so that x becomes the radial variable. The 
solutions of the homogeneous equation are the Legendre 
functions P 3 ^ 2 (/u) and qV 2 {h) and a recurrence relation 
that generates P^ +1 from P^ and P^_ 1 - (Here the order 
fx of the Legendre function has nothing to do with the 
variable /i = yl — kr 2 .) Thus 



Pl /2 (cos X ) = 
Q X J 2 (cos X ) = 



7T 

2 
2 



(sinx) 2 cos 
(sinx)~ 5 sin 



.(4.9) 



+ T(z) (z'riC(z'Mz')dz'. (4.3) P 3 / 2 (cos X ) = 



smx 



v - - ) Pu cos x 



p 



(4.10) 

3/2 

the same relation holds for the Q,/ . It turns out to be 
convenient to write n = v + i. We note that Q2.441 and 
()4.8(l involve this quantity and that n can be real but is 
often imaginary. Thus 



(!) 



sin 3/2 x 



(4.11) 

[cosxcos(nx) — nsinxsin(nx)] , 



similarly writing n when it is real but n — iN when it is 
imaginary: 



Q 3 J-i ( cos x) 



(4.12) 



(^) 2 . V2 [cosxsin(nx) - nsinxcos(nx)] , 
\2/ sin ' x 



3/2 / \ 



— [cos x sinh(iVx) - AT sin x cosh(iVx)] . 

We shall be concerned to have functions which, after mul- 
tiplication by another (sinx) -3 ^ 2 , are nevertheless still fi- 
nite at the origin x = 0. A little expansion around x — 
shows that the P function diverges but Q function satis- 
fies this stringent test. Our next job is to find a solution 
that satisfies this stringent convergence not at x = but 
at the 'other' r = at x = tt- Since that is an alternative 
origin it is clear that [cos(7r— x)] passes that test. A 

little work shows that it is indeed the linear combination 
(2/tt) sin(n7r)P 3 { 2 1 (x) — cos(nir)Q 3 ^\ (x). Finally we no- 

n 2 n 2 

tice that n — 0, which is needed in some of our solutions, 

3/2 

gives Q_ x j2 = 0. This is not a solution at all! However 
lim [(l/n)Q 3 ^ 2 , ,J gives the finite limit 



(2) sin 3 / 2 v [X 



cos x — sin xj 



sin ' x 

We shall therefore use the functions 



(4.13) 



In 



(^dk*^, (4 - 14) 



sm ' x 

Pn = {V) 2 ■ 3/2 S n{K-X) 

\2J sm 3/2 x 

as our independent solutions of the Legendre equation. 
These functions have the added advantage that they re- 
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main real when n = iN 

Sn(x) = 

S'ijv(x) = ™cosx 



sin(nx) 

cosx h smxcos(nx), (4-15) 



n 

sinh(iVx) 
N 



sinxcosh(TVx). 



The Wronskian may be shown to be 
dq n dp n tt sin(n7r) n 2 — 1 



w 



n 1 — /i 2 1 — /j, 2 



(4.16) 



Having formed solutions p and q each of which satisfy 
one of the boundary conditions we look for solutions of 
the inhomogeneous equation of the form 



TD = A{fi)p + B{fi)q. 



(4.17) 



We choose A'p+B'q = 0, and then the equation demands 
that 



(l-^)[A'p , +B , q']=-K, 



(4.18) 



where a dash denotes d/d/i. Solving for A' and B' we 
have, using the Wronskian W/(l — fi 2 ) defined earlier, 
A' = Kq/W. Now p does not satisfy the boundary con- 
ditions at x = 0, so A must be zero there; hence 



.4 



Similarly B' = —Kp/W and to satisfy the boundary con- 
ditions at /i = — 1, x = tt, 

B = -L§ d » = -J"^ m x A *' <420 » 

Thus the solution by variation of the parameters is 



x Kq i""Kp 
P(x) I — sin X d x' + ?(x) I — sin X dx' 



o W 



W 



which gives our solution for uj(x) — (sinx) 3 / 2 cl>: 



(4.21) 



w (x) 



tt/2 



W sin 3 x 



We note that with k = +1 , 4 + TV 2 = A 2 , and 

7r . 9 , sin(n7r) 7T . 9n sinh(A^7r) 
W =2 {n 1)_ n = ~2 (1 + 7V) iV ' ( 4 - 25 ) 



For AT large and x small 
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-—Ne 



-Nit 



1 



A 2 * 1 



(4.26) 

N{it- X )_ 

(4.27) 



Wsin 3 x 3?T'~ V' ' 10 
For N large and x not small nor near tt, 

SiN(x) = ^sinxe Nx , SWO _ x) = ^siu \ < 

Hence our solution near the origin is 

= U 1 + ^f) N f A2fi (x') sin2 X'e-^' d X ', 

(4.28) 

and near the perturbation 

2 Jo N z + 1 Vsmxo/ 

(4.29) 

where at the last line we consider a perturbation with a 
mean fi of H in r ± A with 7VA < 1. 

fc = -1 

The equation to be solved is (|2.43|) with k = — 1 and 
/ = 1. Now we write fx = coshx , iy + \) 2 = A 2 + 4. 
Space is now hyperbolic and fi runs from 1 to oo. The 
relevant solutions of the homogeneous equation are 



p = [Pi 12 + -iQl' 2 



1 /7T 



1 



2V27 sinh 3/2 



-s e (x), 



1 



%Q " 2V2J Sinh 3/2 X 



(4.30) 



where n = (v + 



J5(x) = -(n-l)e- ( " +1)x + (n+l)e- ( "- 1)x , 



S e (x) = -\E{ X ) - E{-x)\ 



Sn(n -X) A 2 fLS„(x') sinx' d X ' The Wronskian 



+ S n ( X ) / A 2 r25 n (7r-x')sinx'dx' 



(4.22) 



dq dp (n 2 — l)n 



^d/Lt d/x 



/i 2 -l 



For x small, 



The solution by variation of parameters is 



1 



sin 3 x 



and for n = iN, 

—7—3 — SiN ~ 

sin x 



(l~n 2 ) 
3 ' 

(1-n 2 ) 
3 



(1 + iV 2 ) 



(l + n 2 )x 21 



10 

(4-n 2 ) 2 
10 X 



(4 + ^V 2 ) 3 
10 X 



1 



n(n 2 — 1) 



p / gAT d/i + q pK d/i 



(4.31) 



(4.32) 



(4.33) 



(4.23) 



(4.24) 



hence, changing the integrations from fj, to x an d cj to w, 
we have 



-3 



(sinhx) 
4(n 2 — l)n 



£(x) / A 2 n(x')^(x') sinhx' d X ' 



(4.34) 



+ S e ( X ) / A 2 n(x , )^(x')sinhx'd X ' 
■>x 
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For small x 

E( X ) = 2-(n 2 - l) X 2 x 



(4.35) 



2n 3n 2 + 1 2 n(n 2 + 1) 3 
~T X+ 12 X 15 X 



so 



Se(x) = f (n 2 -l)x 3 



(n 2 + l) , 
15 X 



. (4.36) 



At large x 

E(x.) = {n + l)e- (n - 1)x = 2(n+l)e- nx smhx, 
S e ( X ) = i(n- l) e ( ,l+1 )* = (n- l)e™ x sinhx.(4.37) 

Near the origin 



1 



(4 - n 2 ) X 2 
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(4.38) 



{n 2 - 4)(n + 1) sinh 2 x'e~" x 'ft(x') d%'. 



At the perturbation 

-to) 4— ff££) e-*'-»'«(x')dx'. 

(4.39) 



2 n J V sinn Xo 



V. THE TIME EVOLUTION OF THE 
DRAGGING 

The evolution of to and Q as functions of cosmic time is 
governed by the equations of motion (contracted Bianchi 
identities) (|2.9|l . For axisymmetric, odd-parity perturba- 
tions these become the angular momentum density con- 
servation law, as discussed in equations 12.22l> - 12.26|) in 
Section 2. In terms of ui(t, r, 0) and Cl(t, r, 9) the conser- 
vation law simply becomes <|2.31ll . i.e. 

[a 5 (p + p){to-tt)]' = (5.1) 
or, in terms of the angular momentum density, we get 

1 j(x,o) 



n — uj = 



a 5 (p + p) r 4 sin 3 ' 



(5.2) 



In this formula the hrst factor singles out the time de- 
pendence of f2 — to. Notice that we have already obtained 
uj(t,r,9) as a function of the angular momentum within 
X, J(< x)i m au three cases k = +1, 0, —1 [see equations 
!pO|) . (|3~14l) ]. We found to to depend on the time 

as l/a 3 (t). Equation (|5.2|l then can be regarded as a 
solution f2(i, r, 6) implied by the equations of motion. 

On the other hand, for Q — lo given at some time t = to 
as a function of Xi^j equation l|5.2() determines the den- 
sity j(X)0) which in turn gives J(< x) an d w(i, x, &) is 
then obtained from equations (|3.4J) . I|3.8|l . ^3.14[l . Angu- 
lar velocity of matter, D(i, x, 0), is then given again by 
equation 15.211 . 

If we are interested in proper azimuthal velocities, we 
can write 



V = ar sin ft. 
and rewrite l|5.2|l as 

V-v 



v = ar sin tt io, 



3{X,6) 



1 (P + P) 



(5.3) 



(5.4) 



Since |fir|, |wr| 1, we have also \V\, \v\ 1. In the 
case of the dust universes (p = 0) the density obeys the 
conservation law pa 3 — constant = C. Equation 1)5. 4|) 
then implies 



V-v 



1 



j(x,Q) 

Cr 3 sin 2 9 a ' 



(5.5) 



This is not valid near t ~ when a — > due to our 
approximation. For a — > oo, V — v — > — the dragging 
becomes perfect. 
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